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Abstract 

We critically review the QCD predictions for the cross sections ai and <jt for diffractive 
p meson electroproduction in longitudinally and transversely polarised states in the HERA 
energy region. We show that both perturbative and non-perturbative approaches which 
involve convolution with the p meson wave function predict values of <jt which fall-off too 
quickly with increasing Q 2 , in comparison with the data. We present a perturbative QCD 
model based on the open production of light qq pairs and parton-hadron duality, which 
describes all features of the data for p electroproduction at high Q 2 and, in particular, 
predicts a satisfactory Q 2 behaviour of (tl/<?t- We find that precise measurements of 
the latter can give valuable information on the Q 2 behaviour of the gluon distribution at 
small x. 
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1. Introduction 



The results of the measurements of p meson electroproduction, 7*p — > pp, are intriguing. 
These are coming from the HI Jl| and ZEUS ||, [| experiments at the HERA electron-proton 
collider, and should be considered in conjunction with the earlier measurements of NMC at 
lower energies. We may briefly summarize the main features of the observed behaviour of the 
cross section cx(7*p — > pp) as follows: 



(i) a ~ 1/Q 5 for 7 < Q 2 < 30 GeV 2 . 

(ii) a ~ W - 8 for 12 < W < 140 GeVQ 

(iii) o L jo T ~ 2 - 4 weakly rising with Q 2 for 6 < Q 2 < 20 GeV 2 . 

(iv) da/dt~e bt with 6 ~ 5 - 6 GeV~ 2 for Q 2 > 10 GeV 2 , 

as compared to 6 ~ 9 GeV -2 for Q 2 = 0. 



As usual, Q 2 is the virtuality of the photon, W is the centre-of-mass energy of the 7*p system 
and t is the square of the four-momentum transfer. The p meson is observed through its 
In decay. If there are sufficient events, then the angular distribution of the decay products 
allows the measurement of the components <jl and ot of the cross section, which describe p 
production in longitudinally and transversely polarised states respectively. As we shall see, the 
measurement of the Q 2 dependence of <jl/<jt is particularly informative. The present data, 
(iii), have large errors, but already indicate the general trend. 

Observations (ii) and (iv) imply the validity of perturbative QCD for the description of high 
energy p electroproduction. Observation (iv) means that the size of the system (the 7* — > p 
Pomeron vertex) decreases with Q 2 , and that at large Q 2 we do indeed have a short-distance 
interaction so that perturbative QCD is justified. In fact the measurement of the slope b ~ 5 — 6 
GeV~ 2 is approximately equal to that expected from the size of the proton, which is consistent 
with the hypothesis that at large Q 2 the size of the 7* — > p vertex is close to zero. From 
observation (ii) we see that the exponent of the a ~ W n behaviour has changed from the 'soft' 
pomeron value n = 4(«p(t) — 1) ~ 0.20 observed in p photoproduction (Q 2 = 0), to a value 
n = 4A ~ 0.8 at high Q 2 which is consistent with the gluon density xg ~ x~ x extracted^ from 
the observed QCD scaling violations of F 2 . Moreover it is in line with the a ~ W ' 8 behaviour 
observed in J/ip photoproduction, where perturbative QCD is expected to be applicable due 
to the sizeable charm quark mass. 

Here we explore the implications of all the observed properties (i)-(iv) for the QCD descrip- 
tion of p electroproduction at HERA. Before we present our detailed study, it is useful to give 
a brief overview of the situation. We begin with the Q 2 dependence of <r(7*p — > pp). We will 

2 This behaviour is observed from the NMC experiment at W ~ 13 GeV right through the HERA energy 
range, 40-140 GeV. 

Corresponding to op(0) ~ 1.08. 

4 The MRS parton sets which best describe the recent HERA measurements of F2 0, and other data, are 
MRS(A') § and MRS(R2) §. For these the effective value of A increases from about 0.2 to 0.3 as Q increases 
from 10 to 50 GeV 2 . 
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show that for p meson electroproduction at high Q 2 , perturbative QCD should be applicable 
to a T , as well as o^. The leading order perturbative QCD prediction for electroproduction in 
longitudinally polarised states is || U 

n(x,q 2 )] 2 m 21 j_ m 

for Q 2 3> m 2 , where x = Q 2 /W 2 and 7 is the anomalous dimension of the gluon density, 
xg(x,Q 2 ) ~ (Q 2 ) 7 . For the relevant range of x, 1CT 3 < x < 10 -2 , we have takenF], for the 
purposes of illustration, the representative average value 7 = 0.3. So the QCD prediction for 
<jl is consistent with the Q 2 behaviour of the data. This is not the case for ot- The prediction 
for ot appears to be too small and to fall too rapidly with increasing Q 2 . If only the leading 
twist component of the light-cone wave functionP] of the p is taken into account then 

m 2 1 

where m is the current (light) quark mass. Although the leading twist is specified by the QCD 
sum rules, the next twist is not known. However, we can make reasonable assumptions to 
estimate its effect. We find that its inclusion has the effect of replacing m 2 in (0) by a factor 
of the order of m 2 . Even considering the uncertainties, the value predicted for <Jl/&t is still 
much too big and has the wrong Q 2 dependence in comparison with the data. We elaborate 
the above arguments in section 2. 

It is frequently claimed that perturbative QCD is not applicable for <jt and that its behaviour 
is of non-perturbative origin, see, for example, ref. 0. But in this case we would expect the 
same slope b as in photoproduction and a 'soft' W ' 2 behaviour. Moreover, non-perturbative 
QCD predicts a 1/Q S or stronger fall-off of ot with increasing Q 2 . Recall that these features 
are not observed in the data. A further discussion of the non-perturbative approach is given in 
section 3. 

Here we present a resolution of the problem, which is based on the application of the hadron- 
parton duality hypothesis to the production of open qq pairs. First we recall the hadron-parton 
duality hypothesis for the process e + e~ — ► hadrons. In this case the hypothesis gives 

- 7* - h)\ ~ (y, < e+e ~ - 7* - > (3) 

/ AM 2 \ 1 / AM 2 

that is the total hadron production (h = p, u . . .) averaged over a mass interval AM 2 (typically 
~ 1 GeV 2 ) is well represented by the partonic cross section. This duality has been checked 




5 From the most recent set of partons [|| [?J we find 7 ~ 0.25 rising to 7 ~ 0.4 as x decreases from 10~ 2 to 
10 -3 for Q 2 w 10 GeV 2 . Of course in the numerical analysis of section 6 the true x and Q 2 dependence of 7 is 
automatically included. 

6 The twist of the p wave function should not be confused with that of the operator which corresponds to 
the 7p amplitude. 
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down to the lowest values of ^/~s. We may therefore expect the duality to apply to diffractive 
p electroproduction for qq produced in the invariant mass interval containing the p meson, 
M 2 < 1 — 1.5 GeV 2 . In this domain the more complicated partonic states (qq + g, qq + 2g, qq + 
qq, . . .) are heavily suppressed, while on the hadronic side the 2n (and to a lesser extent the 
37r) states are known to dominate. Thus for low M 2 we mainly have 

7* -> qq -> 2tt (4) 

or in other words 

( * \ nn V- f M » dcr(7*p^ (qq)p) 2 

a( 7 p - pp) * 0.9 £ 7 m , dM (5) 

g=u,a a 

where the limits M 2 and M% are chosen so that they appropriately embrace the p-meson mass 
region with M 2 — M 2 ~ 1 GeV 2 . The factor 0.9 is included to allow for u production. This 
duality model is predictive. In section 4 we present the QCD formula for open qq electroproduc- 
tion via two gluon exchange, and in section 5 we discuss their general structure. In particular 
we show how the scale dependence of the gluon density softens the ul/o't ~ Q 2 growth with 
increasing Q 2 . The numerical predictions are presented in section 6. There we calculate diffrac- 
tive uu and dd electroproduction and use the duality hypothesis to make detailed predictions 
of the Q 2 dependence of both ai and &t for p meson electroproduction at HERA; results whose 
general structure was anticipated in the discussion of section 5. 

In short, we argue that the convolution of the qq wave function (produced by the 7*) with 
any reasonable p meson wave function would yield a prediction for or which is in disagreement 
with the data. Rather we claim that p electroproduction proceeds via open uu, dd production 
at low M 2 , which has a different structure. Some long time after the interaction with the 
proton, confinement distorts the qq state and forces it to be the p meson, as there are no other 
possibilities. That is the suppression due to the small wave function overlap, (qq\p°), is not 
operative. We depict the situation in Fig. 1. 

2. Standard perturbative approach to the Q 2 behaviour of (Jt,l{p) 

First we wish to sketch the derivation of the perturbative QCD prediction, for <jt shown 
in©, 

<r T (i*p-+pp) ~ [xg(x,Q 2 )] 2 /Q 8 (6) 

for Q 2 ^> m 2 , and to show that it is infrared stable. We must therefore study the 7* — > p 
Pomeron vertex (or so-called impact factor) of Fig. 1(a), which we denote by Jt- We shall also 
consider J L . The factors are given by the convolution of the wave functions if) y (qq) and ip p (qq)- 
It is found that [|, [U] 

Ji = f„J %)Bi (7) 

with % = T or L. The quantity f p is the p meson decay constant and the term e 2 in the quark 
propagator is 

e 2 = z(l - z)Q 2 + m 2 (8) 
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(a) (b) 

Figure 1: Alternative mechanisms for p meson electroproduction: (a) involves a convolution 
of the ip~/{qq) and ip p (qq) wave functions, whereas (b) is based on open qq production and 
parton-hadron duality. At high Q 2 the 'Pomeron' exchange in this picture really stands for the 
exchange of two gluons in the t channel. 

where m is the current quark mass. Bi are the helicity factors}] coming from the quark loop, 
see Fig. 1(a), 

B L = 2z(l-z)yfQ 2 , (9) 
B T = -m. (10) 

if> p (z, kj) is the momentum representation of the p meson wave function; z and are the 
Sudakov and transverse momentum components carried by one of the quarks with respect to 
the photon. The other quark has components 1 — z and — fer- 

The wave functions ip p ,T decrease slowly with k 2 -. and the convergence of the integral in 
(0) is provided only by the denominator e 2 + k\. We therefore introduce an integrated wave 
function 

= f dk 2 T ^ p (z,k 2 T ) (11) 

defined by the scale p 2 = e 2 at which the integral ceases to converge. The quantities 4> p are 
called the leading twist light-cone p meson wave functions and have been well studied in the 
framework of QCD sum rules [12, [T3]. As Q 2 — ► oo (that is e 2 — > oo) we have 



p 



» - 6z(l - z) (12) 



for both i = T, L. Their behaviour at finite scales can be found in Refs. [O, T3|, but in any 
case the <fi p vanish at least as fast clS Z clS Z — ► and 1 — 2! as z — > 1 . We may rewrite the impact 
factors (0) in terms of the integrated wave functions (j> p {z). We obtain 

Ji « /p / f (is) 

7 The vertex satisfies s channel (quark) helicity conservation. In general for f ^ we would also have 
off-diagonal helicity factors, B(^t,Pl) and B(7l,/9t)- 
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Finally, we must convolute J, with the gg-proton interaction amplitude T given by the BFKL 
Pomeron (or two-gluon exchange ladder). The amplitude T behaves as 



hmT = a q -^ p ~ X -^l ~ (e 2 )^ (14) 



where recall that the scale is e 2 = z(l — z)Q 2 + m 2 , and where 7(3;) is the anomalous dimension 
of the gluon. Thus the amplitudes for p electroproduction from transversely (i = T) and 
longitudinally (i = L) polarised photons are 

r dz 

Ai = Ji <g> T = f p J <f> p (z) B t , (15) 

which yield the following Q 2 behaviour of the cross sections 



ax 



\A T \ 2 ~ (m/(Q 2 ) 2 ^) 2 ~ m 2 /Q 6 - 8 , (16) 



ol - \a l \ 2 ~ g 2 (i/(g 2 ) 2 - 7 ) 2 ~ i/g 4 - 8 . (17) 

For illustration, we have again set the gluon anomalous dimension 7 = 0.3. We emphasize that 
the integral in ([TS|) is convergent for At (for any 7 > 0), as well as for Al- Thus e 2 ~ Q 2 
and perturbative QCD is valid not only for ol (where we have additional convergence due to 
B L ~ z(l — z)), but also for a T . 

We note that while the prediction for the relative Q 2 dependence of o~t and <tl is meaningful 
(although not supported by the data), the value for the ratio 

~ L ~ -Q 2 (18) 

(which is in gross disagreement with the data) is not a reliable estimate. The reason is that 
the current u, d quark masses are very small (m < 7 MeV) and that therefore we must consider 
how the non-leading twist contribution to ipj(z, k 2 ^) will modify the prediction for 07. The non- 
leading twist is not known. However, it is reasonable to assume that instead of two variables, 
the p wave function ipj depends on only one variable, namely the invariant mass of the qq pairQ 

where we neglect m 2 . Then, after some algebra, it is possible to show that the impact factor 



Jt can be written in the form of (15) with 0j = 6^(1 — z), and that the helicity factor becomes 



B T = -\ m p [z 2 + (1 - zf] , (20) 

rather than the very small 'leading- twist' prediction given in (|10|) . The reason that we still 
obtain a definite prediction for Jt, again in terms of f p , is due to the fact that this same non- 
leading twist component of ipj describes the decay px — > e + e~, that is the kx integral over the 



3 This hypothesis is very natural from a dispersion relation viewpoint ma 
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quark loop describing the px decay is the same integral that occurs in the impact factor J? for 
Q 2 ^> m 2 . In this way we are able to normalise the non-leading twist to the observed width of 
the decay, that is to the decay constant f p . 

If we estimate the p electroproduction amplitude At of ( |T5| ) using the modified form ( |20| ) 
of Bt then we obtain 

o T m 2 

(21) 



°L " Q 2 



with c ~ 2. The precise value of c depends on the actual forms of (p^ ,L (z) at the exper 



p 



imentally relevant scales, p 2 ~ 10 GeV 2 , which are far from the asymptotic region where 
4^' L ( Z ) — 6^(1 — z). In our approximate estimate of c ~ 2 we have used the 4>^' L (z) wave 
functions of ref. |13| . Although a considerable improvement on (|18D, the prediction (^1|) for 
the ratio a T /a L is still much smaller than the observed ratio, and as before decreases more 
rapidly with Q 2 than indicated by the data |], ^ f|. In short, the standard perturbative QCD 
predictions for <Jt{i*P —>■ pp) are not in agreement with the observations. 



3. Non-perturbative approach to the Q 2 dependence of (Jt{p) 



It has been argued that the main contribution to Of comes from the non-perturbative re- 
gion p. Let us disregard the fact that the perturbative integral ([15]) is convergent for a T an d 
suppose that non-perturbative effects dominate. In order to obtain non-perturbative contribu- 
tions associated with small (~ p 2 ) virtualities, we must get contributions from the end-point 
regions of integration 



< 



Only then will we sample small scales e 1 



and 1 — z < p 2 /Q 2 . 
p 2 and large distances p r 



(22) 

1/e ~ 1/yU. However, for 



large distances the quark effectively has a constituent mass m q ~ ^m p and the non-relativistic 
wave function, 0j(z) ~ S(z — |), is appropriate. Certainly 4>^(z) decreases exponentially, or at 

1. Thus the contribution from the regions (22) should 
v z{l — z), as in (|T2"D, we would obtain from ( |T3| ) with 



least as a large power, as z —>■ or z 
be strongly suppressed. Even if 4>^(z 



e ~ p 



a T (non-pert. 



p 2 /q 2 



dz 



z)Bn 



1 



(23) 



Thus for the actual non-perturbative prediction we would expect an even faster fall-off with 
increasing Q 2 . 



4. QCD model for (Jl,t{p) via open qq production 

The above discussion suggests that the problem in successfully describing p meson electro- 
production may be associated with having to convolute with a p meson wave function, which 
inevitably leads to a form-factor-like suppression of the form \ (qq\p°}\ ~ 1/Q 4 - Here we study 
an alternative and physically compelling mechanism for p electroproduction based on the pro- 
duction of uu and dd pairs in a broad mass interval containing the p meson. In this mass 



6 



interval phase space forces these qq pairs to hadronize dominantly into 2ir states, with only a 
small amount of 3tt production. Moreover, provided the gg-proton interaction does not distort 
the spin, we expect that the process 7* — > qq — > 2ir will dominantly produce 2tt systems with 
J p = 1~. The calculation of the diffractive electroproduction of qq pairs therefore allows, 
via the parton-hadron duality hypothesis, a detailed prediction of the structure of p meson 
electroproduction. 




Figure 2: Diffractive open qq production in high energy 7*p collisions, where z is the fraction of 
the energy of the photon that is carried by the quark. The transverse momenta of the outgoing 
quarks are ±A/r, and those of the exchanged gluons are ±£y. 



The formula for the diffractive production of open qq pairs is given in ref . |T5|, [Ljj . For light 
quarks we may safely put the current quark mass m — 0. The process is shown in Fig. 2. We 



use the same notation as in ref. HT5| , so the scale at which the gluon distribution is sampled is 
denoted 

K 2 = z(l - z)Q 2 + kl 



k T 



where the last equality follows since z(l — z) 

P = 



k 2 T jM 2 and 

Q 2 



Q 2 + M 2 ' 



(24) 



(25) 



Note that the scale K 2 plays the role that e 2 played for exclusive vector meson production (cf. 
(|H)), and that it determines the transverse distances bx ~ 1/K that are typically sampled in the 
process. It is convenient to replace the dk\ integration over the quark transverse momenta k? 
in formulae (40, 41) in ref. [T5[] by an integration over dK 2 . Then it is straightforward to show 
that these formulae giving the JltP ^ (qq)p cross sections in the forward direction (t = 0) 
may be written in the form 



d 2 cr L 
dM 2 dt 

d 2 a T 
dM 2 dt 



Air 2 e 2 a 



Q 2 



An 2 e 2 q a 



(Q 2 + M 2 ) 4 
M 2 



dK 2 K 2 



K 2 



3 (Q 2 + M 2 y 



r(Q 2 +M 2 ) 



K 2 



1 - 4K 2 /(Q 2 + M 2 ) 
dK 2 (l - 2f3K 2 /Q 2 ) 
1 - AK 2 /(Q 2 + M 2 ) 



[h(K 2 )Y, (26) 



[I T {K 2 )Y (27) 



where a is the electromagnetic coupling. The quantities I^t are the integrations over the 
transverse momenta, ±£t, of the exchanged gluons (see Fig. 2) 





= I<2 J 


di\ 


h(K 2 ) 




r dt% 



a s {e T ) f{x,t T ) — - — I (28) 



± 




\K 2 


Kj) 


(± 


1 


\K 2 




) 2 " 





K 2 - 9k 2 4- 



where x = (Q 2 + M 2 )/W 2 , 

K 2 = ^J(K 2 + t 2 T ) 2 - Ak 2 £ 2 T , (30) 

and f(x,£^) is the unintegrated gluon distribution of the proton. We will use formulae ( |26| ) 
and ( p?D to predict p meson electroproduction. They involve integration over the quark fcf (or 
i^ 2 ) and over the i\ of the exchanged gluons. As we are dealing with a diffractive process we 
see that the cross sections have a quadratic sensitivity to the gluon density. 

It is useful to inspect the leading lni^ 2 approximation to the dGp integrations of (|28|) and 
(|29|). In this approximation it is assumed that the main contributions to the integrals come 
from the domain &p < K 2 , and so on expanding the integrands we obtain 

/- = = s^p f f f(x , e T) = ^1 xg{x , *•>. (3i) 



By analogy with flT4]), we see that Il,t are essentially the cross sections for the (qq)L,T in- 
teraction with the proton. Of course, in the calculations presented in section 6 we do not 
use the leading log approximation, but instead we perform the explicit d&p integrations over 
f(x,£j,) = d(xg(x, 4))/<91n£y given in (p8|) and (j29|). We treat the infrared region using the 
linear approximation described in ref. [Hi] for low t\ values (that is l 2 T < Q. We find stability 



of the results to reasonable variations of the choice of 



5. Insight into the structure of the cross sections gl,t 

In section 6 we show the predictions for the Q 2 behaviour of <jl and ot for p electroproduc- 
tion, which are obtained from the numerical evaluation of ( p6|) and ( p7|) integrated over the p 
mass region. However, it is informative to anticipate some of the general features of the results. 
First we study the infrared convergence of the dK 2 integrations of ( p6|) and (|27|). We note from 
the approximate forms of Il,t in fl3T|) that 

h ~ x- x {K 2 y/K 2 (32) 

where A and 7 are the effective exponents of the gluon defined by 

xg(x,K 2 ) ~ x~ x (K 2 y. (33) 

We see that the integration (|26D is infrared convergent provided that 7 > as K 2 —> 0, whereas 
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K 2 = 2.5 GeV 2 
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Figure 3: The continuous curves show the effective anomalous dimension 7 of the gluon (defined 
by xg(x, K 2 ) ~ (if 2 ) 7 ) determined from the MRS(R2) set of partons for K 2 = 2.5, 5, 10 and 
20 GeV 2 . The dashed curves correspond to the values of 7 for the R4 set of partons. 



we require 7 > 0.5 to ensure the convergence of (|27|). How does the value of 7 depend on K 2 1 
At high energy W (that is x ~ Q 2 /W 2 — ► 0) the gluon g(x, K 2 ) increases much faster as x 



decreases for large K 2 (xg 



~ x 



with A > 0.3) than for small K 2 . Thus the effective anomalous 



dimension 7 increases when x and K 2 decrease. The behaviour is evident in Fig. 3 which shows 
the values of 7 (as a function of x for selected K 2 ) obtained from two recent sets of partons. 
For example, let us take a typical value x ~ Q 2 /W 2 = 10~ 3 relevant for the measurements at 
HERA (say Q 2 = 10 GeV 2 and W = 100 GeV). We see from Fig. 3 that 7 increases from 0.3, 
0.45, 0.6 to 1 as K 2 decreases from about 20, 10, 5 to 2.5 GeV 2 . The infrared convergence 
requirement, 7 > 0.5, of ( |2T| ) is therefore already satisfied when K 2 has decreased to 8 GeV 2 . 
In general, the behaviour of 7 with K 2 amply provides, via (0), the infrared convergence of 
(|2"7D , as well as of (^). This explains the reason why our numerical evaluation of (E7|) for ot 
depends only weakly on the infrared cut-off^ Kq. Indeed integral ( p6|) for <tl is controlled by 
contributions close to the upper limit and we expect 



d 2 a L 



(Q 



2\27-2A 



dM 2 dt Q 6 (34) 

This is exactly the same Q 2 behaviour as the prediction (0) for exclusive pl electroproduction; 
9 In section 6 we choose Kq = 0.2 GeV, the order of the inverse confinement radius 1 fm _1 . 
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here we have been more precise and displayed the Q 2 dependence coming from the x(~ Q 2 /W 2 ) 
behaviour of the gluon. Of course, the result ( p4[ ) is very approximate and the detailed depen- 
dence of the Q 2 behaviour of <jl (as well as ax) on the properties of the gluon must await the 
numerical predictions of section 6. 

Nevertheless we can take the general discussion further and anticipate the main features of 
the Q 2 behaviour of the important ratio <J L /ax- We first rewrite (^) and (^) in terms of an 
integration over the angles of the produced qq pair. We use the polar angle 9 of the outgoing 
q in qq rest frame with respect to the incident direction of the proton. Thus we have 



h 



and the square root in the denominators of 
the numerator of (WK) 



|Msin0, 



(35) 



]) and ( p7|) is equal to cos 9. Also the factor in 



1 - 2(3K 2 /Q 2 = Ul + cos 2 



\d\M\ 2 + 



where d Xfl (9) are the conventional spin rotation matrices. Then equations (p6|) and 



d 2 a L 
dM 2 dt 

d 2 a T 
dM 2 dt 



47r 2 e 2 a 



Q 2 



47r 2 e 2 a 



(Q 2 + M 2 ) 2 
M 2 



dcos9 



-i 



3 (Q 2 + M 2 ) 2 4 



dcos9 



-i 



d\ (9) 



d\i(e) + dU{9) 



(36) 
become 

(37) 
(38) 



where the dependence on the rotation matrices appropriately reflects the decay of the p meson 
from longitudinally and transversely polarised states, respectively. 

In the limit of no interaction with the proton (that is II = It = constant) the photon has 
to produce the qq pair in a pure spin J = 1 state. We immediately find from fl3~T|) and ([!8|) that 

a L Q 2 J d cos 9 sin 2 9 Q 2 

&T 



(39) 



2M 2 Jdcos9(l + cos 2 9) AM 2 ' 

In the realistic situation, the two-gluon exchange interaction distorts the qq state produced by 
the 'heavy' photon. Some idea of the consequences of this distortion can be anticipated from 
the leading log approximation ( j3l|) for II and It, in which 

(K 2 V 1 



K 2 



sin 



1 1-7 ' 



(40) 



We substitute this behaviour into ( |3"TD and 
the underlying qq production amplitudes (~ 
We then use the identity 

' * sin p 9 d9 



, and project^ out the spin 1 components of 



d\ x (9) 1(9) where I L 

nh + b) 

r(i + \ P ) 



1(9) 



sin 



(41) 



10 To be precise the rotation matrices D J 



—(f)) form the orthogonal basis and we project out the com- 
ponents c(A) from the qq amplitudes D\* X I{9) with the matrix D\ x . However, the integrations are trivial and 
hence the projection can be done simply in terms of d\ x . 
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to evaluate the projections 

2J+ 1 



c(A) = — ±_ J dcos6 (d{f(6)I(6)) d{f(9), (42) 

assuming that 7 is a constant over the region of integration. With this assumption we find the 
interesting result 



&l ____ Q 2 |c(A = 0)p 



a T 2M 2 |c(A = 1)| 2 + |c(A = -l)p 




(43) 



The dependence on 7 has the effect of masking the Q 2 growth of <Jl/&t- This can be seen by 
inspecting Fig. 3 - higher Q 2 means larger x and both changes imply smaller 7. The projection 



integrals fl42 ) for the amplitudes (with their linear dependence on Ii(6)) are more infrared 
convergent than fl3~T| ) and (|38|). Now ot (as well as <Jl) is convergent provided only that 7 > 
as K 2 — > (that is as 6 —>■ 0). In fact, provided x remains sufficiently small, both the <7l and 
<jt integrations receive their main contributions from the region K 2 < Q 2 /4, and so we should 
insert into (f43"D the average 7 sampled in this x, K 2 domain. Indeed the decrease of 7 with 
increasing K 2 < Q 2 /4 is found to considerably suppress the growth of ctl/o't with increasing 
Q 2 , and to largely remove the gross disagreement of the QCD prediction with the data; see the 
full numerical calculation presented in section 6. We may turn the argument the other way 
round. Accurate measurements of the ratio ctl/^t a s a function of Q 2 will offer an excellent 
way of constraining the K 2 and x behaviour of the gluon g(x,K 2 ) in the region K 2 < Q 2 /4 
and x ~ Q 2 /W 2 . Of course result fl4*3|), which is based on a constant 7, is oversimplified. It is 
given only to indicate the general trend. The full calculation of section 6 is performed with a 
realistic gluon distribution and so automatically allows for the K 2 (and x) dependence of 7. 

We see that the projection integrals (^) converge in the infrared region of small K 2 « 
^Q 2 sin 2 9 (that is at small 9) for any 7 > 0, even for ot (that is for c(A = ±1)). We have 
stronger infrared convergence for gl or c(A = 0) due to d\ Q = — sin#/\/2. We also notice that 
the factor 1(9) = l/(sin 2 #) 1-7 , arising from the gg-proton interaction, gives a strong peak in 
the forward direct ionF]. It means that the distortion caused by the interaction will, in principle, 
produce higher spin qq states. Most probably the higher spin states at small M 2 are killed by 
confinement during the hadronization stage as there is insufficient phase space to create 2ir 
states with large spin with M 2 < 1 GeV 2 . In any case the higher spin components^ cannot 
affect p production, since confinement cannot change the spin of the produced qq state. At 
higher energies (small x) the anomalous dimension 7 grows and the function 1(8) is not so 
singular as 9 — > 0. Therefore in this energy domain the incoming spin of the qq system is not so 
contaminated by J 7^ 1 components arising from the interaction with the proton. In the black 
disk limit of the proton when the cross section approaches the saturation (unitarity) limit 7 
tends to 1 and we come back to pure J = 1 qq production. 

n The height of the peak is limited by the infrared cut-off, Kq = 0.2 GeV, provided by confinement. 

12 Indeed it will be interesting to study the detailed spin decomposition of 7* — > open qq production as 
a function of M 2 . In this way we can investigate how the QCD 'Pomeron' distorts the initial state and how 
confinement/parton-hadron duality operates in different (relatively small) M 2 regions for the different J p states. 
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6. Numerical QCD predictions for p electroproduction 

We use parton-hadron duality to predict p electroproduction from the QCD formulae for 
open uu and dd production. To be precise we compute 

/•(1.05 GeV) 2 (j(j T rr( J —\\ 

°lAp) = 0.9 / ( dM 2 daL ' T Z~ 1} (44) 

J (0.6 GeV) 2 dM 2 

where d<7L,T(J = l)/dM 2 are the spin 1 projections of open qq production of fl3~T|) and (p8|), 
carried out as described in fl4"2"|), and where the cross sections have been integrated over t 



assuming the form exp(— b\t\) with the observed slope b = 5.5 GeV -2 [|l], ||. The factor 0.9 is 
included in (Q) to allow for u production. The II and It integrations over the gluon transverse 
momentum are computed from (|28| ) and (|29|) as described in ref . [ ll| . We checked the stability 



of the results to contributions from the infrared regions of the dK 2 and d£ T integrations. First 
we varied the infrared cut-off around the value K Q = 200 MeV that we used to evaluate ( f2~6|) 
and fl2~T|). Second we explored the effect of varying i 2 , around the value = 1.5 GeV 2 that we 
used to evaluate the integrals of (|28|) and (|29|). Recall that we use the linear approximation 
described in ref. J0| to evaluate the contribution from the region £ T < t\. We found only a 
weak sensitivity to variation of the choice of £q. For instance reducing i\ to 1 GeV 2 changes 
the cross sections by less than 5%. We will report on the sensitivity to variation of K at the 
end of the section. 

We begin by taking the gluon distribution from the MRS(R2) set of partons 0, which 
correspond to a QCD coupling which satisfies as(Mf) = 0.12. The parton set with this QCD 
coupling, found by global analysis of deep inelastic and related data (including recent HERA 
measurements of F 2 ), is favoured by the Fermilab jet data with E? < 200 GeV 0. We first 
compare our cross section predictions obtained with this gluon with the data. Then we use 
different gluon distributions from several recent sets of partons to study the sensitivity of 
7*p — > pp to the behaviour of the gluon. 

Note that we use phenomenological gluon distributions which are obtained from global 
fits to deep inelastic experimental data, rather than "ab initio" distributions calculated from 
theoretical models. Thus the gluon distributions that we use already incorporate absorptive 
effects. 

There is another crucial ingredient in the calculation of the cross section for diffractive 
open qq production. Virtual gluon corrections to the process shown in Fig. 2 are surprisingly 
important. The relevant diagrams are discussed in ref. [JTSJ and lead to tt 2 enhancements of 
the 0(as) corrections. If the contributions are resummed they lead to an enhancement of the 
lowest order result by a factor exp(asCV7r), the so-called K factor enhancement, where the 
colour factor Cp = 4/3. A similar K factor is well known in Drell-Yan production, although 



there the contributions come from different virtual diagrams |L5[ . For the Drell-Yan process the 
enhancement can be as much as about a factor of 3. In our case the K factor can, at present, 
only be estimated. It proves to be the main uncertainty in the normalization of diffractive qq 
production. The major ambiguity is associated with the choice of the argument of as- We 
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• HI: 40 GeV < W < 140 GeV 
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Figure 4: The predicted Q 2 dependence of the cross section for 7*p — » pp compared with (a) 
HI data H collected over the energy range 40 < W < 140 GeV and (b) preliminary ZEUS data 
in energy bins with (W) = 56, 81 and 110 GeV. The QCD curves for the various values of 
W are obtained using MRS(R2) partons M. 



take the scale to be 2K 2 . Since the K 2 integrations are dominated by contributions towards 
the upper limit this choice is equivalent to a scale < Q 2 /2. With this choice we obtain the 
values of the 7*p — > pp cross section shown by the curves in Fig. 4, which are in reasonable 
agreement with the measured values. For our choice of scale the average K factor for <jl varies 
from about 3 to 3.7 for Q 2 going from 25 to 10 GeV 2 , and is about 20-25% larger for <jt (as in 
this case somewhat lower K 2 values are sampled). The cross section agreement shown in Fig. 4 
corresponds to a physically reasonable choice of scale, and leads to a sensible range of size of the 
K factors. It shows that the open qq duality model for p electroproduction is at least consistent 
with observations. Due to the sensitivity to the choice of scale, clearly the agreement cannot 
be regarded as confirmation of the approach. Nevertheless, it does imply the existence of a 
sizeable 'V 2 " enhancement of the Born amplitude, as was also found in the Drell-Yan process. 

On the other hand the predictions for the Q 2 dependence of the ratio <Jl/&t have much 
less ambiguity. The calculations are compared with the measurements at HERA in Fig. 5. 
The agreement with the data shows a dramatic improvement over the QCD expectations which 
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R = c L /c T 
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Figure 5: The Q 2 dependence of the QCD predictions for the ratio u^jaT of the electropro- 
duction of p mesons (7*p — > pp) in longitudinal and transverse polarisation states compared 
with the most recent HI and ZEUS || data. MRS(R2) partons are used. 

involve convolution with the p meson wave function. The small x behaviour of the gluon plays 
a crucial role in masking the Q 2 increase anticipated in these earlier predictions of the ratio. 

The dependence on the gluon is seen in Fig. 6 which compares the Q 2 behaviour for o~l/o"t 
at W = 90 GeV for the gluon distribution of several recent sets of partons (MRS (A') ||, GRV 
JT7[ , MRS(R2) 10]). We stress that the normalization of the QCD predictions for the cross 
section are dependent on the choice of the mass interval embracing the p meson and on the 
estimate of the K factor enhancement. On the other hand the ratio <Jl/&t is not so sensitive 
to these ambiguities. At this stage it is relevant to study the stability of the results to variation 
of the infrared cut-off K . This we also show in Fig. 6, where we present QCD predictions 
based on MRS(R2) partons for two different choices of Kq. We see that the cross section is 
hardly changed while the ratio <Jl/&t increases a little when K is increased from 200 to 300 
MeV. Such a result is to be anticipated as ot samples, on average, smaller K 2 values than 
<7l- However, we see that the sensitivity of the predictions for a^jor to the value of K is 
sufficiently weak so that measurements of the ratio can give a reliable probe of the gluon. 
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a(y*p->pp) [nb] R = o L /o T 




Figure 6: The QCD predictions for W = 90 GeV based on three recent sets of partons P, |7|, |I7|1 
compared with the recent HERA data JT|, [|. We also show the sensitivity of the predictions 
using the MRS(R2) partons to the choice of the cut-off K Q ; the dot-dashed curves correspond 
to Kq = 300 MeV whereas all other curves correspond to Kq = 200 MeV. The dot-dashed curve 
in (a) essentially coincides with the continuous curve which demonstrates the insensitivity of 
the cross section prediction to the value of K , whereas we see that the ratio cr L /a T of (b) has 
some dependence. 



7. Conclusions 

We have shown that the diffractive electroproduction of p mesons at high Q 2 can be described 
by perturbative QCD. Indeed, since p production in both longitudinally and transversely po- 
larised states is being measured at HERA with better and better precision, the process 7*p — > pp 
can serve as an excellent testing ground for QCD. Moreover, we have shown that it also provides 
a sensitive probe of the small x behaviour of the gluon distribution. 

The validity of perturbative QCD is ensured by the large value of Q 2 . This is already 
suggested by several features of the existing data [Q, ||, ||. However, the measurements of the 
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ratio cr L /a T do not support the behaviour, 



OL Q 2 



(45) 



(Jt 2m 2 ' 



predicted from QCD by convoluting 7* — > qq diffractive production with our knowledge of the 
p meson wave function. The main problem is that the predictions for <jt are too small and fall 
off too quickly with increasing Q 2 . We showed that a non-perturbative approach to <jt does not 
resolve the conflict with the data. Rather we argued that on account of the low mass of the p 
meson the convolution with the wave function should be omitted. The uu or dd pairs produced 
in the p mass region have, because of phase space restrictions, little alternative but to hadronize 
as 27r states. Thus a more appropriate approach to p electroproduction is to apply the parton- 
hadron duality hypothesis to open uu and dd production. Indeed we found that this model gives 
a good description of all of the features observed for diffractive p electroproduction at HERA, 
including in particular the Q 2 behaviour of o^jaT- To gain insight into expectations of the 
model, we first made a simple estimate based on assuming a constant anomalous dimension 7. 
We found 



where M is the invariant mass of the qq pair and 7 is the effective anomalous dimension of the 
gluon defined by xg(x, K 2 ) ~ (K 2 ) 1 , where the typical K 2 sampled is K 2 < Q 2 /4 (approximated 
to be the same for both and aj). The decrease of 7 with increasing Q 2 masks the strong 



idea of the crucial role played by the gluon distribution. In Figs. 4-6 we showed the results 
of the full calculation. The computation is based on a measured gluon distribution and so 
automatically allows for the appropriate K 2 and x dependence of 7. The figures compare the 
detailed predictions of the model with the measurements of diffractive p electroproduction at 
HERA. The main uncertainty is in the normalization of the cross section. One source is in the 
choice of the width of the AM 2 interval over which to apply the duality hypothesis. The second 
is associated with the K factor enhancement which arises from virtual gluon corrections to open 
qq production. The normalization is sensitive to the choice of scale used as the argument of 
as in the calculation of the K factor. The data show evidence for a K factor of about 3-4, 
comparable in size to the K factor enhancement established for Drell-Yan production. 

The QCD model prediction of the ratio ctl/ct is essentially free of the above ambiguities. 
Fig. 6 shows that precise measurements of the ratio for p electroproduction at different values 
of Q 2 , and the 7*p cm. energy W, will provide a valuable probe of the behaviour of the gluon 
distribution g(x, K 2 ) in the kinematic domain x ~ Q 2 /W 2 and K 2 < <5 2 /4. 




(46) 



growth shown in (pE5|). Of course result (SET) is greatly oversimplified but it gives a good 
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